Electronic-vibrational coupling in single-molecule devices 
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Experiments studying vibrational eflFects on electronic transport through single molecules have 
observed several seemingly inconsistent behaviors, ranging from up to 30 harmonics of a vibrational 
frequency in one experiment, to an absence of higher-harmonic peaks in another. We study the 
different manifestations of electronic-vibrational coupling in inelastic and elastic electron transport 
through single molecules. For the case of inelastic transport, higher harmonics are shown to be 
damped by additional small factors beyond powers of the electron-vibration coupling constant A. 
Two mechanisms greatly increase the size of secondary peaks in inelastic transport: coupling between 
electron transport and spatial motion of the molecule, and the "pumping" of higher vibrational 
modes of the molecule when vibrational excitations do not completely relax between electron transits. 



Electronic devices made from single molecules offer 
both a new means of observing fundamental physics, such 
as the Kondo eflFect 0, , and the potential to circum- 
vent the device size limitations of silicon electronics 0. 
Transport in molecular electronic devices is strongly af- 
fected by several mechanisms that are negligible in ordi- 
nary semiconductor electronics. An example is the ob- 
servation of molecular vibrational modes in DC trans- 
port measurements on single-molecule devices made from 
fuUerenes Q and other organic molecules [^]. Coupling 
between molecular vibrations and the electronic states 
responsible for transport was also proposed as an expla- 
nation for the strong temperature dependence of the cur- 
rent through single molecules 0. 

There are also experiments, however, in which the 
electronic-vibrational coupling appears quite weak, as 
in the STM measurements on an adsorbed molecule re- 
ported in 1^. This paper discusses the most important 
mechanisms for electronic-vibrational effects in transport 
through single molecules. The fact that different exper- 
iments probe different mechanisms explains much of the 
observed variation. The apparent coupling measured in 
an experiment depends on whether elastic or inelastic 
tunneling is being probed; whether transport electrons 
"pump" the molecule into higher vibrational states; and 
whether the molecule's motion affects its coupling to the 
leads. Our calculations are based on single-level Einstein- 
type models for electron-phonon coupling, and in what 
follows both "phonon" and "vibrational mode" describe 
the vibrational states of the active molecule. 

The chief simplifying assumptions in our calculations 
are that the system is at zero temperature (T — 0) and off 
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resonance, so that the conductance per molecule g <^ 
For the systems under consideration the phonon energy 
huj and bare level width T are both typically larger than 
the experimental temperature. Under the assumption 
that incident electrons reach the molecule in its ground 



state, the vibrational modes of the molecule appear first 
as inelastic steps in the differential conductance (when 
an electron passing through leaves behind one or more 
phonons), and then as subpeaks of resonant-tunneling 
(elastic) peaks. The height of inelastic steps is shown 
to fall off much more rapidly with phonon number n 
than the weight of elastic subpeaks, explaining the lack 
of higher-harmonic peaks in |3] . This results from an or- 
thogonality constraint and hence should hold beyond the 
simple Einstein model used here. 

The transport of electrons becomes more complex 
when the phonons excited by one electron's passing 
through do not have time to relax before the next electron 
arrives: there is a nonthermal steady-state of phonon 
occupancies on the molecule, and high phonon-number 
states are "pumped" as soon as the bias voltage is larger 
than the phonon energy. We also consider a model where 
electron hopping on/off the molecule is modified by vibra- 
tional excitations, motivated by the experiments 
The final part of the paper discusses which predictions 
of the idealized Einstein-type models can be expected to 
apply more generally. Quantities such as overlap inte- 
grals and electron-phonon coupling strength, which here 
are taken as inputs, can be calculated ab initio with mod- 
erate accuracy Q for specific molecular systems. 

The Einstein Hamiltonian for an electronic level of en- 
ergy Erf and a single vibrational mode of energy Tilij is 



Xnujc''c{b + b''). 



(1) 



This Hamiltonian conserves electron number: the zero- 
electron state with p phonons, |0,p), has energy phco. 
The one-electron ground state can be found by a simple 
canonical transformation and has energy id = ed — 
X'^hiv. It is annihilated by 6-1- A and is a superposition of 
states of different phonon number: 



(0,n|c|l,0) 



A"e-^'/2 



(2) 
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Straightforward algebra gives the matrix elements be- 
tween general zero- and one-particle eigenstates: 
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(0,n|c|l,fc) = (n|^(A + 6^)'=e-^ /2 V-^^|m) 



k-l 



/-f m— n 

e-'^^ E L. (3) 
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To describe an idealized single-molecule device, the 
hopping of electrons on and off the molecule must be 
added: 



^{tk2cliCk2 + h.c). 



fcl 



(4) 
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Here i?Li, are the Hamiltonians of the left and right 
leads, and tki is the overlap integral between state fcl and 
the molecular level. This hopping model is appropriate 
when the molecular level is weakly coupled, but the low 
conductance of current devices suggests that this is a 
valid assumption It is assumed below that the leads 
are essentially uniform in energy over the ranges probed 
in an experiment, and that the lead densities of states 
Pl.r and hopping elements ' are energy-independent. 
The left and right level widths are 'yL,R = '2ttpl,r\V^'^\'^ . 

Steps in differential conductance appear in the Ein- 
stein model when the bias voltage is a multiple of the 
phonon energy. When the bias voltage V > nhu, an elec- 
tron passing through the molecule can excite n phonons. 
We find that the step height falls off very rapidly with n 
for reasonable parameters because of an orthogonality re- 
quirement, which may explain the absence of n > 1 peaks 
in inelastic STM spectroscopy of single molecules, j] We 
start by assuming that between electron transits, any 
phonons which may have been excited by the first elec- 
tron decay without electron transfer before the second 
electron arrives. Later this assumption will be weakened, 
with dramatic effects on the step strength. In the limit of 
rapid phonon decay, each electron reaches the molecule 
in its ground state |0, 0). The current for id > 2V occurs 
via virtual excitation of the molecule: 



'n<eV/huj 



n=0 

2erirfl 



7i<eV/huj peV/2 



2-Kh 



E 



A{n, E) 



E 



„_g J —eV/2+nhu) 

(0,n|c|l,fc)(l,fc|ct|0,0)' 
ed + khu — E 



A{n,E) dE, 



(5) 



When eV crosses nfioj, the contribution of the new chan- 
nel /q^k is initially zero. / is continuous and dl/dV 
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FIG. 1: Inelastic steps in differential conductance: (dl/dV) 
in units of Go = 2e'^TLTR/{2-Kh{nijf), versus V{hLo)-^, for 
X = 2, id ~ 3.75hu!. Even at this strong coupling, higher steps 
are much smaller than the first step because of additional 
factors of [hu/ed). 



has steps at eV = nhu for zero temperature, which 
broaden with increasing temperature but remain visible 
for kT < fiLo. The magnitude of these steps can be found 
explicitly: define the dimensionless step height 



h 

2^ 



dl 
dV 



\eV—nhuj~'^ 



di_ 

Iv 



\eV—nhuj 



(6) 



and the function f{x) = {-X^ [Ti-x) - T{-x, -A)] 
E^-a^TT^. Then 



11=0 n\{x+n) ' 

Ki = 

Ko = 



r,r 



2ir(nui) 



[f{ld/nu - i)- 
f{idlfiu + \)f 



^X^e--^-\fi~ed/nc.- 

2f{id/nuj) + fiid/huj + 1) 



(7) 



and the pattern of discrete derivatives continues to higher 
steps. The first two step heights are plotted in Fig. 1, 
normalized by the smooth part of the differential conduc- 
tance at eV = fiuj. 



G= 4h 



dV \ eV=huj 



2e^ TlT 
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(8) 



The main qualitative prediction of the above is a sim- 
ple consequence of the orthogonality of phonon states: 
inelastic steps are weak (compared to inelastic peaks) 
because the sum in © vanishes when {Tiui/ed) — > 0. For 
small A and small fku/ed, the height of the nth step de- 
creases with n as {Xhcu / id)"^" ■ Even for A = 2 the n > 1 
steps are very small, as seen in Fig. ^ The elastic tun- 
neling peaks in the Einstein model (discussed below) fall 
off only as A"; the additional factor hui/id is absent. 

The subsidiary steps can be much larger if the phonon 
decay time is slow compared to the current, as then 
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phonons are predominantly excited and relaxed by the 
transport electrons. Before considering this "pumping" 
mechanism, we briefly review the effect of phonons on 
elastic peaks. Under the weak-coupling assumption of 
this paper, the differential conductance is simply propor- 
tional to the electron density of states on the molecule. 
Start by fixing A = so that there is one resonant tun- 
nehng peak, at eV/2 = id- The electron density of states 
is Lorentzian, and Ap per spin integrates to one electron: 



Ap{E) 



The differential conductance is then 



dl 
dV 



2e2 



Li R 



ApieV/2) 



n (ri + r^) 



(9) 



(10) 



and the total peak weight is Wq = 2eTi,Tii/h{TL + Tr). 

The above results generalize quite simply to the Ein- 
stein model as long as the level broadening , Tr <C 
huj. The nth resonant peak is a Lorentzian centered on 
eV/2 = id + nhcu and has width (F^ -I- rii)/2, as before, 
and weight W = Wq I (0, 0|c|l, n)p under the stan- 
dard resonant tunneling assumptions for an interacting 
level [nl in [13 ■ Hence the total weight in all the dif- 
ferential conductance peaks is the same as before; the 
broadened lineshape gene rated by electron-phonon cou- 
pling is discussed in [lj|. These elastic subpeaks of a 
resonant tunneling peak have been observed in 0. For 
small A, the peaks fall off as A^"; note that the factor 
{hio/id)'^"' which was present for the step heights in the 
inelastic case is now absent. 

The current through a single-molecule device is 
strongly modified when vibrational excitations on the 
molecule do not relax completely between electron tran- 
sits. There is a crossover between two regimes with in- 
creasing current: low current is described by the preced- 
ing results, while at high current the vibrational excita- 
tions are predominantly excited and relaxed by transport 
electrons, rather than by decaying directly into the leads. 
Calculating the crossover requires knowing the rate of vi- 
brational excitations from decay into the bulk without 
electron transport, which depends sensitively on the de- 
tails of the molecule-lead contacts. The phonon lifetime 
for nanoparticles on surfaces is discussed in [lij |. 

We consider for brevity the limit where vibrational 
excitations on the molecule are created and annihilated 
only by transport electrons. This situation is expected 
to be more easily realized using molecules on surfaces [^ ; 
molecules strongly contacted to leads will have rapid de- 
cay of phonons into the bulk leads. The rate equations 
obtained in this limit can be easily modified to include vi- 
brational relaxation without electrons; there is a smooth 
crossover from decay entirely through transport electrons 
to decay entirely through phonon relaxation into the bulk 



leads. The matrix element for a typical process in second- 
order perturbation theory is 



fiEL)^tLtRj2 



(0,j|c|l,n)(l,n|ct|0,fc > 
ed — X^hoj — El + {n — j)fujj 



(11) 

This notation means that the electron moves from L to 
R and starts with energy El, in the process changing the 
vibrational state from j to k. 

The j ^ k terms vanish if u}/(ed — El) — > 0, since 
then the denominator is constant and the sum evalu- 
ates to 5jk/{id — El)- The pumping mechanism depends 
on having a vibrational mode frequency not too much 
smaller than the energy scale of the electronically ex- 
cited state. On the other hand, a large vibrational fre- 
quency means that peaks only appear at higher voltages, 
so there is a window of vibrational frequencies which 
will have the most significant effects. This may partly 
explain the experimental observation '3| that only one 
or a few vibrational modes seem to show up strongly 



in transport. As required by unitarity. 



fiEL 



^R^L 



{EL + {k- j)nuj). 



The total rate for the molecule to change from state j 
to state fc, fc > j, is then 



T^L^R. 



dEL TLTn 



(0,j|c|l,n)(l,n|ct|0,A;> 



ed 



X^hu) — El + {n — j)fiLO 



.(12) 



Here Tl,r = 2TT\tL,R\^ Pl,r and pl,r is the density of 
states (assumed constant as a function of energy). In 
what follows symmetric coupling is assumed: F = F^ = 
Ffl. Similar equations with more terms apply for the re- 
laxation rates, since there are more ways the molecule can 
relax. These rate equations are similar to those for spin 
excitations in [l^ but with many more states involved. 

The rate equations for virtual (second-order) tunnel- 
ing cease to apply once sequential (first-order) tunneling 
becomes energetically allowed. Sequential tunneling ef- 
fectively closes the rate equations, since high vibrational 
levels of the molecule are quickly relaxed once sequential 
tunneling becomes possible [nfioj -f /i^ > Id)- This is 
again a consequence of our assumption that the molecule 
is weakly contacted, i.e., F ^ T'^ /ed- As an example of 
the pumping effect, we calculate that for A = 1, = 5.5, 
the first step is enhanced by a factor 2.49 relative to no 
pumping, and the second step by a factor 7.23. These 
numbers include no transport by sequential tunneling 
from highly excited states; this would further increase 
the step size in the pumped case. Now we return to 
rapid relaxation (no pumping) but consider how spatial 
motion of the molecule modifies tunneling. 

It has been suggested 0, [^ that the vibrational modes 
observed in experiments are related to the center-of-mass 
oscillation of the molecule. Such molecular motion affects 
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the overlap of the electron wavefunction on the molecule 
and the leads, and hence the sequential tunneling rates 
and Ffl. To study its effect on the peaks in differen- 
tial conductance, we use a simple generalization of the 
Einstein model, where Tl and F/j are now dependent on 
the oscillator state n. Assume that the the electronic 
state is a Gaussian centred on the molecule, and that 
the electronic states on the leads are i/^l = ^/2/roe~^^^" , 
V'fl — 1/27^06'-''"^-'/'''' where the leads arc at a; = and 
d, and tq is the spread of the wavefunction. 

The wavefunction of the molecule in the nth harmonic 
oscillator state is 
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where I — y^Ti/muj, Hn are Hermite polynomials and 
xq is the mean position of the molecule. In a Born- 
Oppenheimer-type approximation, discussed below, the 
overlap integrals are obtained and averaged over the 
probability distribution of the molecule in the nth har- 
monic oscillator state. The resulting tunneling rates are 

= (V^)e'-^Ao%^Vn^e-^-/-L„(^) (14) 
ro 

where L„ are Laguerre polynomials and is the spread 
of the electronic state. Note that for a molecule in a very 
tight potential, i.e. I — *■ 0, the results are independent 
of n. The initial increase in Ti^^jj with n results in an 
increase in the peaks observed in the differential conduc- 
tance, and depends only on the relative spread of the 
wave functions at the leads and that of the harmonic os- 
cillator states. Fig. □ shows the resulting current-voltage 
curves with and without spatial dependence for moder- 
ate electron-phonon coupling A = 1. The higher har- 
monic peaks have enhanced weight and show more than 
seven appreciable steps; these results agree well with the 
increase in the strength of higher peaks observed in ^ . 

The Born-Oppenheimer approximation that gives H14fl 
is justified in the high-current limit, where a large number 
of electrons tunnel through the molecule in a single pe- 
riod of the oscillator. In the case of very low current and 
temperature (for example, less than one electron in every 
100 periods tunnels through the molecule in j^), the anal- 
ysis is more complicated because it is no longer justified 
to average over the molecule's probability density. The 
finite temperature and molecule-environment interaction 
in real experiments will help validate Born-Oppenheimer 
approximation of incoherent molecular motion. 

The above results explain, first, why inelastic tunnel- 
ing experiments often fail to see higher phonon peaks, 
and second, why elastic tunneling experiments see a large 
number of peaks even for moderate electron-phonon cou- 
pling. This paper constructed an explicit model of direct 



FIG. 2: Steps in current-voltage due to elastic tunneling pro- 
cess (a) without and (b) with dependence of tunneling rates 
r^jPfl on vibrational state of the molecule. In (b), the pa- 
rameters are A = 1 and l/ro = 1; in (a), there is no spatial 
dependence {I 0) and again A = 1. 



coupling between molecular motion and electronic tun- 
neling, which quantitatively explains the current-voltage 
characteristics seen in experiments 0, with reason- 
able values of the electron-phonon coupling and vibra- 
tional amplitude. Strong enhancement of vibrational ef- 
fects also occurs when the current is sufficiently high for 
"pumping" of phonons by transport electrons. 
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